Exact ground state properties of the one-dimensional Coulomb gas 
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The ground state properties of a single-component one-dimensional Coulomb gas are investigated. 
We use Bose-Fermi mapping for the ground state wave function which permits to solve the Fermi 
sign problem in the following respects (i) the nodal surface is known, permitting exact calculations 
(ii) evaluation of determinants is avoided, reducing the numerical complexity to that of a bosonic 
system, thus allowing simulation of a large number of fermions. Due to the mapping the energy 
and local properties in one-dimensional Coulomb systems are exactly the same for Bose-Einstein 
and Fermi-Dirac statistics. The exact ground state energy has been calculated in homogeneous and 
trapped geometries by using the diffusion Monte Carlo method. We show that in the low-density 
Wigner crystal limit an elementary low-lying excitation is a plasmon, which is to be contrasted with 
the large-density ideal Fermi gas/Tonks-Girardeau limit, where low lying excitations are phonons. 
Exact density profiles are confronted to the ones calculated within the local density approximation 
which predicts a change from a semicircular to inverted parabolic shape of the density profile as the 
value of the charge is increased. 

PACS numbers: 71.10.Pm Fermions in reduced dimensions (anyons, composite fermions, Luttinger liquid, 
etc.) - 71.10.Hf Non-Fcrmi-liquid ground states, electron phase diagrams and phase transitions in model 
systems - 73.21.Hb Quantum wires 



The recent progress in nanoscale technology made it 
possible to realize clean one-dimensional quantum gases, 
such as ultracold atoms confined in elongated traps 
electrons in single- well carbon nanotubcs|3 and in semi- 
conductor quantum wires Q. A peculiarity of a onc- 
dimensional world is that such systems can not be ex- 
plained by the conventional Landau theory of normal 
Fermi liquids, but rather by an effective low-energy Lut- 
tinger liquid description^, generalizable to long-range 
Coulomb interactions 0. A number of approaches of in- 
creasing accuracy (random-phase approximation, Singwi- 
Tosi-Land-Sjolander scheme, density functional theory) 
have been devised to study energetic and structural 
properties 0. The most precise calculations of energy 
have been obtained by the Monte Carlo technique as in 
Ref. where a quasi-one-dimensional geometry with a 
finite width of the transverse confinement is studied. The 
purpose of the present work is to carry out exact calcu- 
lations in a strictly one-dimensional geometry. 

In this paper we use the Bose-Fermi mapping to find 
the ground state energy of a one-dimensional single com- 
ponent Coulomb system exactly within statistical pre- 
cision. The mapping applies both to homogeneous and 
trapped geometries. 

We consider a single-component system of TV particles 
(bosons or fermions) of charge e and mass m in a one- 
dimensional box of length L. Periodic boundary condi- 
tions are applied and the Coulomb potential is truncated 
for interpartcle distances larger than L/2. The Hamilto- 
nian reads as 



Natural length scales in a homogeneous system are de- 
fined by atomic units, that is Bohr radius uq = jme^ 
for length and the Rydberg Ry = jla^ for energy. The 
system properties are governed by a single parameter, 
fixed by the ratio of Bohr oq and Wigncr-Scitz = l/2n 
radiuses or, equivalcntly, by the dimensionlcss density 
riflo with n = N/L. 

Hamiltonian ([1]) describes a system with the pure 
Coulomb interaction VintiT") = e^/r between charges not 
screened while the movement of charges confined to a 
quantum wire, so that excitations of the levels of the 
transverse direction is not possible and the system be- 
haves as quasi-one dimensional. For a harmonic trans- 
verse confinement this implies that the energy associated 
with the Hamiltonian ([1]) is small compared to NH'^/mb'^, 
where b is the transverse oscillator length. The pecu- 
liarity of the Coulomb interaction in a one-dimensional 
geometry is that its Fourier transform cannot be calcu- 
lated due to a divergence at short distances. Commonly, 
this ultraviolet divergence is cured by considering the 
ground state of the transverse harmonic oscillator. This 
effectively introduces a cutoff to the interaction poten- 
tial and removes the divergence as V^J/ (0) = y/ne'^/b. 
Here, instead we note that the unbound potential in ^ 
will not produce any divergence in the energy while the 
pair distribution function vanishes when two particles ap- 
proach each other. This permits us to successfully apply 
the one-dimensional Bose-Fermi mapping @ and relate 
the many-body ground-state wave functions of fermions 
ipp{zi, • • • , zn) and bosons ^'^(^i, • • • , zn) as 
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In the case of bosons the divergent 1 /r interaction mimics 
the Pauli exclusion principle, which for fermions prohibits 



2 



any two fermions from overlapping. 

We resort to the diffusion Monte Carlo (DMC) tech- 
nique for finding the ground state properties of the sys- 
tem. The guiding wave function is chosen in bosonic Bijl- 
Jastrow form, adding antisymmetrization in the case of 
fermions: 



V'_f(zi, 



<zn) = H-^^d^i " Zj|)sgn(zj - Zj) . (3) 



It is well-known that simulations of fcrmionic systems 
suffer from the "sign problem". In terms of the DMC 
algorithm this means that a "fixed-node" approximation 
should be used, that is the fermionic ground state wave 
function is expressed as some non-negative function mul- 
tiplied by a function which has a certain nodal surface. 
The DMC method finds the best energy for a given choice 
of the nodal surface, providing an upper bound for the 
ground state energy. The exact ground state energy is 
reached only if the nodal surface is exact. Importantly, 
the nodal surface of Eq. ([3]) is exact. This means that 
the mapping ([3]) permits to find the exact energy of the 
fermionic system in DMC calculation. Furthermore, the 
same mapping allows avoiding evaluation of determinants 
and reduces the complexity of the calculation from 
to N'^, making accessible a larger number of particles. 
We choose the two-body correlation term f2{z) as 
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where coefficients Ci , C2 are fixed by the continuity con- 
ditions for f2iz) and f2iz) at the matching distance Rpar- 
The short-range part of f2{z) corresponds to the zero- 
energy solution of the two-body scattering problem. This 
avoids divergences in the energy as two particles come 
close to each other. The long-range part of f2{z) is taken 
from a hydrodynamic theory The parameter Rpar is 
optimized by minimization of a variational energy. 

The typical kinetic energy per particle scales quadrat- 
ically with the density, opposite to the linear density de- 
pendence in the Coulomb interaction energy. At large 
densities the kinetic energy becomes dominant and can 
be approximated by the energy of an ideal Fermi gas 
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In the opposite regime of low density the energy can be 
obtained by summation of the Coulomb potential energy 
under the assumption that all particles are localized at 
lattice sites of a Wigner crystal. The leading contribution 
to the energy (sec, for example, Ref. [Toj ) 
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for a fixed number of particles is linear in the density. In- 
stead, for a fixed density the energy per particle and the 
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FIG. 1. (Color online) Energy per particle as a function of 
the dimensionless density for A*' = 1000; 100; 10 (from top to 
bottom) particles. Circles, Monte Carlo results; dash-dotted 
line IFG/TG energy, Eq. ((5|; dashed lines, energy of a Wigner 
crystal, Eq. ((SJ. 



chemical potential diverge logarithmically as the num- 
ber of particles increases. This is a consequence of the 
long-range nature of the Coulomb interaction potential. 
In a three-dimensional system a similar divergence can 
be cured by the "jellium" model, that is, by restoring 
charge neutrality by adding a uniform charge of the op- 
posite sign. In a one-dimensional system introduction 
of a uniform opposite charge would lead to a logarith- 
mic divergence in the short-range part of the Coulomb 
interaction potential. 

The energy of a homogeneous system is shown in 
Fig. [TJ At low densities the energetic properties are 
precisely captured by the Wigner crystal description. 
At the same time no true crystal is formed as quan- 
tum fluctuations destroy the diagonal long-range order 
in a one-dimensional geometry [llj. The large-density 
regime is well described by an ideal Fermi gas (IFG) or 
Tonks-Girardeau (TG) gas depending on particle statis- 
tics. There is a smooth crossover between the Wigner 
crystal and IFG/TG gas regimes. 

The low-lying excitations in gases with short-range 
potentials are phonons with a linear dispersion rela- 
tion ujk = |fc|c, where c is the speed of sound. This 
is the case of (5- interacting Lieb-Liniger [l^l, Calogero- 
SutherlandfTsj. and Tonks- Girardeau |Q gases, in which 
the long-range properties of the correlation functions are 
described by the Luttinger liquid theory Q. The long- 
range nature of the Coulomb interaction potential leads 
to different low-lying excitations which instead are plas- 
mons. Their dispersion relation can be found from the 
classical equations of motions in a Wigner crystal[i,[l3. 
The low-momentum dependence is nonanalytic and con- 
tains logarithmic terms 



2e^n 



In^ 

n 



H 



(7) 



When the excitation spectrum is exhausted by one 
branch of excitations, the excitation spectrum is related 
to the static structure factor Sk = {p\pk) /N by the Feyn- 
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FIG. 2. (Color online) Static structure factor for = 100 
particles in the Wigner crystal regime (nao — 0.01). Symbols, 
DMC results; lines, Feynman relation Sk ~ h'^k^/{2mujk) 
for different dispersion relations: phonons and plasmons, 
Eq. (0; harmonic theory [hw{k)/Ryf = 4^(3) -2ij3(e^'=/") - 
2Li3(e-"=/") H. 
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FIG. 3. (Color online) Density profile of A'' = 10 trapped par- 
ticles. Solid line, diffusion Monte Carlo result; local density 
approximation with equation of state as in ideal Fermi gas 
(dashed line) and Wigner crystal (dash-dotted line). DMC 
density profiles are calculated using the technique of pure es- 
timators [l^ . 



man formula Sk ~ fi^k'^ /{2mujk)- Wc perform a DMC 
calculation of the static structure factor and present the 
results in Fig. [5J We observe that the low-energy charge 
excitations are well described by the plasmonic dispersion 
relation ([7]). We note that although the logarithmic term 
in ([7]) mathematically is the leading one, it alone is not 
enough to describe accurately the excitation spectrum 
and the use of the subsequent term is very important. 
In the opposite regime of large density the dispersion re- 
lation is similar to that of an ideal Fermi gas, and the 
low- lying excitations are phonons. 

The long-range nature of the Coulomb potential causes 
no divergences in finite-size trapped systems. We con- 
sider the case of a harmonic oscillator and add the con- 
finement term Vtrapiz) = to Hamiltonian ([TJ. 
It is convenient to use oscillator units of hio and oscillator 
length ttosc = y^h/mtd for energy and distances, respec- 
tively. In dimensionless units the Hamiltonian reads as 
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and the system is characterized by the number of parti- 
cles N and the dimensionless charge q = \J aosc/a,Q- 

For a small number of particles the density profile has 
a typical shell structure of a mesoscopic system as shown 
in Fig. [3l When the shell structure can be neglected, 
the system properties are well described within the local 
density approximation (LDA) The chemical poten- 
tial /i of the trapped system is then approximated by the 
sum of the external potential Vtrap{z) and the chemical 
potential ^hom{n) of the homogeneous system, with the 
local density n{z). The value of the chemical potential is 
fixed by the normalization condition N = J n{z) dz. 

In the limit of small charge g — > 0, the chemical 
potential of a homogeneous system as in Eq. can 
be approximated by a kinetic energy of an ideal Fermi 
gas ^hom = 7r^?i^n^/2m and is independent of particle 



charges. The density profile is a semicircle of a Thomas- 
Fermi size R = V2Naosc'- 



(9) 



The chemical potential is /i = Nhiu and the energy per 
particle is E/N = Nhjj/2. Even if the LDA does not 
describe the shell structure, thus leaving out contribu- 
tions of the order of 1/iV, the obtained LDA result for 
the energy turns out to be exact for a trapped ideal Fermi 
gas. 

In the limit of large charge q — ^ oo, the chemical po- 
tential of a homogeneous system has a linear dependence 
on the density ^hom = ^e^nlnN , which is obtained by 
differentiating the Wigner crystal energy ({6]). The den- 
sity profile has an inverted parabola shape, typical for 
gases with short-range interaction potentials |15|. 
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with Thomas-Fermi radius R = (3iV In 7Vq2)i/3_ xhe 
chemical potential is /x = 1/2 (Sg^iVlniV)^/^ and the 
energy per particle E/N = 3/10 {3q^NlnNy/^. 

We adapt the guiding wave function ([3]) of a homo- 
geneous system by multiplying it by one-body terms 
fi{zi) = exp(— aZj^) with the free parameter a opti- 
mized by minimizing the variational energy. The two- 
body terms are chosen as /2(z) = ^/zli{2,jqz). 

Figure U] shows the energy of N charged particles in a 
harmonic trap. For small charges the energy per par- 
ticle is well described by that of an ideal Fermi gas 
and increases linearly with the number of particles N . 
Each new particle occupies the next I = N — 1 level 
of the harmonic oscillator and the energy increases by 
El = {I + l/2)huj. For large charges the local density ap- 
proximation result with Wigner crystal equation of state 
applies. In this regime the energy per particle increases as 
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FIG. 4. (Color online) Ground-state energy of iV = 100; 10; 5 
(from top to bottom) charged particles in a harmonic trap. 
Circles, DMC results; dashed line, energy if IFG/TG gas; 
dash-dotted line, LDA result for Wigner crystal. 



(A^ln A^)'^/^. In between the two limits there is a smooth 
crossover. 

The density profile for a system of TV = 10 is presented 
in Fig. [3] For small charges the LDA predicts a semicircu- 
lar shape of a radius which depends only on the number 
of particles and is independent of the value of the charge. 
For that reason different density profiles with ^ 0.1 
almost coincide. The regime of large charges is a trap 
analogue of a Wigner crystal. The particles are highly 
localized and form well defined shells, while the overall 
shape is an inverted parabola of a width which depends 
on the value of the charge. Differently from a homoge- 
neous Wigner crystal, here the distance a between lattice 
sites is no longer constant and depends on the distance z 
from the center of the trap as a{z) w l/n(z). The LDA 
description is known to be imprecise at the borders of the 
system [m as it predicts cusps in the density, while the 
true density profile is smooth as can be seen in Fig. [31 
Furthermore, the large-density Wigner crystal equation 
of state is not applicable at the edges where the density 
is small, which worsens the LDA description of a trapped 
Wigner crystal at the edges. The square of the size of the 



system is related to the potential energy of the trap and 
has a similar behavior to the total energy, see Fig. 

To conclude, we have studied properties of a single 
component Coulomb 1/r gas in a one-dimensional ge- 
ometry. The short-range divergence of the interaction 
potential permits us to exploit the Fermi-Bose mapping 
for the ground state wave function and to solve the Fermi 
sign problem in the considered system as (i) the mapping 
allows exact evaluation of the ground state properties, (ii) 
the numerical complexity of a simulation in a fermionic 
system is the same as in a bosonic system. This per- 
mits us to carry out ah initio simulation of many-body 
systems with as many as 1000 fermions. Due to the map- 
ping, the bosonic and fermionic Coulomb systems have 
the same energy and local properties (density profile, pair 
distribution functions, etc.). 

We use the diffusion Monte Carlo method to do an ex- 
act calculation of the ground state energy of a system in a 
box with periodic boundary conditions and in a harmonic 
trap. In the regime of large density a homogeneous sys- 
tem is well described by an ideal Fermi gas for fermions 
or Tonks- Girardeau gas for bosons. There is a smooth 
crossover transition to the Wigner crystal regime as the 
density is lowered. We show that the low-lying excita- 
tions are plasmons as manifested in the low-momentum 
behavior of the static structure factor. 

The energy and the density profile in trapped sys- 
tems are obtained within the local density approxima- 
tion and are confronted with the exact results of the 
diffusion Monte Carlo calculation. In the ideal Fermi 
gas/Tonks-Girardeau limit the system properties are de- 
fined by the number of particles and are independent of 
the charge. The limit of large charge is a trapped ana- 
logue of a Wigner crystal where localized shells play a 
similar role to crystal lattice sites. The amplitude of the 
oscillations in the density profile experiences a drastic 
change in the shape varying from a semicircle (IFG/TG 
limit) to an inverted parabola (Wigner crystal limit). 
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